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We analyze the strong coupling limit of spin-one bosons in low dimensional lattices. In one- 
dimensional lattices, for an odd number of bosons per site (No = 2k + 1), the ground state is a 
dimerized valence bond crystal with a two-fold degeneracy; the low lying elementary spin excita- 
tions carry spin-one. For an even number of bosons per site (iVo = 2k), the ground state is a 
nondegenerate spin singlet Mott state. We also argue that in square lattices in a quantum spin 
disordered limit, ground states should be dimerized valence bond crystals for an odd integer No. 
Finally, we briefly report results on exotic states for non-integer numbers of bosons per site in one- 
dimensional lattices. 



PACS number: 05.30.Jp, 05.50.+q, 75.10.Jm, 75.45.+j 



Recently, Greiner et al. reported the observation of 
superfluid and Mott insulating states in optical lattices 
[1]. The experiment convencingly demonstrates the fea- 
sibility of investigating many-body states of cold atoms 
by varying laser intensities of optical lattices [2,3]. Be- 
cause of controlable properties of optical lattices, it is 
now widely accepted that many strong-coupling limits 
which are not accessible in solid-state systems can be in- 
vestigated in optical lattices. Furthermore, since optical 
lattices are free from imperfections, an accurate compar- 
ison between theories and future experiments appears to 
be likely. Consequently, we believe that studies of spin- 
one atoms in optical environments might also create new 
opportunities for understanding some fundamental issues 
in quantum magnetism. This in turn should be useful to 
possible applications toward quantum information stor- 
age and processing. 

Spin correlated condensates have been recently investi- 
gated in experiments [4,5]; two-body scatterings between 
spin-one atoms were found to result in either ferromag- 
netic or antiferromagnetic condensates [6-8]. Due to 
spin-phase separation there also exist a variety of novel 
condensates which carry distinct fractionalized topolog- 
ical excitations [9-11]. Interesting spin correlated frac- 
tional quantum hall states were also analyzed [12,13]. 

For spin-one bosons with antiferromagnetic interac- 
tions in high dimensional optical lattices, spin correlated 
Mott insulating states were investigated very recently 
[10]. Spin singlet and nematic Mott insulators were 
pointed out for interacting spin-one bosons in strong- 
coupling and intermediate coupling limit respectively. 
In this Letter we examine spin correlations in low di- 
mensional Mott insulating states. We show that in the 
strong-coupling limit in one-dimensional optical lattices, 
there are two distinct Mott states depending on the 
" even-odd" parity of numbers of spin-one bosons per site 
{No). When Nq is an even number, the ground state 
is a non-degenerate spin singlet Mott state and the low- 
est lying elementary spin excitations correspond to spin 
5 = 2 ones; for an odd integer No, the ground state is a 
dimerized valence bond crystal which has a two-fold de- 



generacy and an elementary spin excitation in this state 
corresponds to a spin-one state. We also report results 
for one-dimensional lattices with non-integer numbers of 
atoms per site and square lattices. To get oriented, let us 
first consider the following coupled "mesoscopic" spinor 
condensates. 

Two coupled condensates of spin-one atoms 

Consider atoms such as sodium ones which interact with 
each other via antiferromagnetic interactions [6-8]. As- 
suming all atoms condense to a single one-particle or- 
bital state, we have shown that the dynamics of this con- 
densate is characterized by a constrained quantum rotor 
model [9-11]. The effective Hamiltonian of two coupled 
spinor condensates (at site A and B) of spin-one bosons 
with antiferromagnetic interactions can therefore be ex- 
pressed as 
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In Eq.l, at each site k two unit vectors nfc (defined on 
a two-sphere) and e %Xh (defined on a unit circle) char- 
acterize the orientation of an 0(3)- and 0(2)- rotor re- 
spectively; they are conjugate variables of the spin oper- 
ator Sfe and atom number operator Nk, i.e. [Sfecn^^] = 
—iheapyTLik-y, [Nk,Xk] = ih- The total spin of a conden- 
sate at site k is therefore the angular momentum of the 
0(3)-rotor at the same site. 

TLhop originates from hopping of bosons from one site 
to other and conserves the total spin of two coupled 
condensates, t is equal to Noto] to is the hopping 
integral of an individual boson, t can be varied by 
changing laser intensities in optical lattices [1]. E s = 
{2IY 1 = ATrh 2 {a 2 - a )p /3MN and and E c = (2C) _1 
= 47rfi, 2 (2a2 + ao)po/3MNo characterize the spin and 
"charge" gaps in the excitation spectrum. They depend 
on the atomic mass M, the atomic number density po 
and two-body scattering lengths 02.0 in the total spin 
S = 2,0 channels for two spin-one cold atoms. And 
for sodium atoms, a 2 exceeds ao by a few percent [6]. 
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N n (^> 1) is the number of bosons per site and /j is the 
chemical potential. 

Finally, the Hilbert space of the Hamiltonian in Eq.l 
is furthermore subject to a constraint (— l) Sfc +" fc = 1 so 
that the corresponding microscopic wave functions are 
symmetric under the interchange of two bosons; Sk,rik 
are the eigenvalues of Sk,Nk respectively (see section II 
B in Ref. [11] for more discussions on the constraint). 
We will consider a situation where E s is much smaller 
than E c because the difference between a 2 and ao is 
much smaller than ao as for sodium atoms; furthermore 
in the strong- coupling limit which interests us, t is much 
smaller than both E s and E c . 

When the energy gap E c (^> E s ) is much larger than 
the hopping energy t, bosons are localized at each site; 
and if there are No = 2k bosons at each site and the hop- 
ping is absent (t = 0), the ground state for each isolated 
site is an 5 — state where each pair of atoms forms 
a spin singlet, following Eq.l and the constraint. This 
limit was also discussed in [8]. Excitations have to be 
spin S = 2,4, 6... states and are gapped by an energy of 
order of E s . Therefore, two weakly coupled condensates 
also have a spin singlet ground state and all excitations 
are gapped by either E s or E c . Particularly, there will 
be no excitations at an energy scale much lower than E s . 
In lattices of all dimensions, ground states in this limit 
are spin singlet Mott insulating states as pointed out in 
a previous work [10] (see Fig. la)). 

The situation for an odd number of particles per site is 
more involved. When the hopping is absent, the ground 
state for each individual site is a spin-one state with a 
three- fold degeneracy. Eq.l and the constraint further 
indicate that excitations are states with spin S = 3,5, ... 
which are gapped from the spin-one ground state. So at 
energies smaller than E s , each individual spinor conden- 
sate can be treated as a spin-one pseudo-particle. For two 
decoupled sites, the ground state has a nine-fold degen- 
eracy, corresponding to states with total spin 5 = 0, 1, 2. 

A detailed calculation, taking into account virtual hop- 
pings between two sites, shows that the effective ex- 
change interaction (for E s smaller than E c ) is TLab = 
- J Es=o,i.2 a sPs(S A + S B ). Here P S (S A + S B ) is the 
projection operator of a total spin 5 state of the two 
condensates A and B. J = t 2 /E c is the strength of the 
exchange energy and is much smaller than E CjS ; 
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where 7l (0) = 1/6, 73 (0) = 2/15; 7s (l) = 1/10; 72 (2) = 
2/15, 73 (2) = 7/150; 72 (0) = 7 i(l) = 72 (1) = 7 i(2) = 0. 
c s = E s /E c ; and ao is always larger than ol\^- 

Eq.2 indicates that two coupled spinor condensates, 
each of which has A^o = 2k + 1 atoms again have a spin 
singlet ground state. The lowest lying spin excitations 
correspond to 5 = 2 ones, whose energies in this case 



are equal to e(S = 2) = j(«o — a 2 ). In addition, there 
are excitations of spin 5=1, with energies e(5 = 1) 
= J{a - «i) and e(5 = 1) > e(5 = 2). In the follow- 
ing we will focus on ground states and excitations in low 
dimensional optical lattices with odd numbers of bosons 
per site in the strong-coupling limit. The "even-odd" 
feature emphasised in this letter is absent in the phase 
diagram for spinless bosons discussed in Ref. [2] . 

ID lattices with odd numbers of bosons per site 
For the study of one-dimensional lattices, we generalize 
the two-site Hamiltonian in Eq.l to 



H 1D = - J ^2 asPsiSi + Sj). 



(3) 
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The sum is over < ij >, all pairs of nearest neigh- 
boring condensates in onc-dimcnsional lattices. Eq.3 is 
equivalent to the bilinear-biquadratic spin chain model 
H/J = cos?jX;<i,> s » • S J + sin? ?E<y>( S i • s i) 2 ; if V 



satisfies tan r\ 



{ol\ 



c\2 — 4ao)/((*i — a 2) and the sign 



of sin 77 is chosen to be the same as that of a\ — ^a 2 — gao- 
The ground state of the bilinear-biquadratic Hamil- 
tonian depends on r\ as suggested in [14]. The 77 = 
point represents a usual 5 = 1 Heisenberg antiferromag- 
nctic spin chain which was also studied in [15]. When 
77 varies between — 37r/4 and — 7r/4, one expects that the 
ground state of coupled condensates should be a dimer- 
ized valence-bond crystal, or DVBC. 

As c s defined after Eq.2 varies from to +00, 77 varies 
from 7/0 to — 7r/2; 770 = — 7r/2 — arctanl/2. For spin-one 
bosons such as sodium atoms, c s is much smaller than 
unity and our calculation shows that r\ (« 7/0) is precisely 
within the dimerized regime. This is further supported 
by a direct mapping between the problem of interacting 
spin-one bosons and a quantum dimcr model (see be- 
low and Ref. 16 for more discussions); it is also consistent 
with destructive interferences between Z 2 instantons in 
one-dimensional lattices for an odd number of atoms per 
site [16]. A straightforward calculation yields that the en- 
ergy per site of a DVBC state is f = -[§(a + ^) + ^], 
with ao,i,2 given in Eq.2. 

In a dimerized configuration shown in Fig. lb), each 
two neighboring condensates are in a spin singlet state, 
or form a valence bond; this results in a valence bond 
crystal doubling the period of the underlying optical lat- 
tices. A DVBC state can be considered as a projected 
BCS pairing state of pseudo-particles introduced above. 
For one-dimensional lattices of 2M (M — > +00) sites with 
a periodic boundary condition, a DVBC state reads (up 
to a normalization factor) 

\9 >= Vg.gV 2 m cxp ( 2 g a * + (Qo = an)) |0 > (4) 
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Aq m as a bosonic creation operator of a spin-one pseudo- 
particle with a crystal momentum Q and m = 0, ±1; 
h(Q) = exp(iQ) and Q G [— 7r,7r]. go = 1, ,9i = ±1 cor- 
respond to two-fold degenerate DVBC states [17]. Vg.g. 
is a generalized Gutzwiller projection to forbid a double 
or higher occupancy of pseudo-particles at each site and 
Vim is a projection into 2M-particle states. Eq.4 repre- 
sents a projected superposition of singlet pairing states in 
(Q, —Q) and (Q, ir — Q) channels. It appears to be plau- 
sible to observe interference between these two channels 
[18]. 
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Fig.l Schematic of microscopic wave functions for differ- 
ent one-dimensional Mott insulating states of spin-one 
atoms, a) A spin singlet Mott state for an even num- 
ber of atoms per site; b) a DVBC Mott state for an 
odd number of atoms per site; as a reference, we also 
show in c) a nematic Mott insulating state which is sta- 
ble only in high dimensional lattices [10]; d) a spin-one 
excitation (kink-like) in a DVBC state; the location of 
the kink is indicated by a green dashed line. In a), 
b), and d), each pair of blue and red dots represents 
a spin singlet state of two spin-one atoms. In c), d), 
a unpaired blue dot carrying an arrow pointing toward 
the direction of ft represents an atom in a spin-one 
state specified by fi = (0,0): ^= sin0exp(— i(f>)\l, 1 > 
— sin#exp(i(/>)|l, — 1 > + cos#|l,0 >; in c), localized 
atoms at each site "condense" in an identical spin-one 
state. 



The DVBC state breaks the crystal translational sym- 
metry and has a twofold degeneracy; one of the degen- 
erate states is characterized by the following correlation 
functions (i ^ j) 



<S i .S i >=--(l + (-l)'> i)i+lj 

< S* • S i+1 Sj ■ S j+1 >=i(l-(-l) i )(l-(-l) i ). (5) 

Eq.5 indicates long-range order in valence bonds; the al- 
ternating feature between site i — 2k + 1 and i = 2k in 
the correlation function reflects a two-fold degeneracy of 
DVBC states [19]. 

An S = l,S z = rn excitation can be created by the 
following product operator 



?ec T 



(6) 



wh e re *+ = ^(A+A+_ x + A%_ X A% X - A+ A+ ) is a 
creation operator of a valence bond at link £ which con- 
nects two neighboring sites i and j; A+ m is the creation 
operator of a pseudo-particle or a collective state at site 
7 with spin S = 1, S z = m. Finally, the product is car- 
ried over all links £ along path C 7 which starts at site 7 
and ends at the infinity; furthermore, C 7 is chosen such 
that the first link along the path is occupied by a va- 
lence bond. Eq.6 also represents a spin-one domain wall 
soliton(of S z = 0, ±1) located at site 7 in the DVBC as 
shown in Fig. Id). Taking into account the hopping ma- 
trix of the domain wall %j = — ^cto[5i,j-2 + ^j.j+2] along 
one-dimensional lattices, we obtain the following band 
structure for spin-one excitations: 
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with Q x defined as a crystal momentum of the excitation, 
— 7r/2 < Q x < 7r/2. These excitations are purely mag- 
netic and involve no extra atoms, i.e. they are "neutral". 

ID lattices with non-integer numbers of bosons 
per site When the number of bosons per site devi- 
ates from an integer, extra bosons can hop along one- 
dimensional lattices and form a superfluid state. The 
situation with Nq = 2k ± e, e <C 1 is particularly simple. 
In this case, 2k spin-one bosons localized at each site are 
in a spin singlet state; extra bosons added to or removed 
from the system tend to condense in a state with zero 
crystal momentum. However, hyperfine spin dependent 
scatterings destroy spin long-range order of condensed 
spin-one bosons in one dimension as pointed out in [9,11]. 
The ground state therefore is spin disordered but phase 
coherent, similar to that in one-dimensional traps. Inter- 
actions between bosons furthermore lead to usual Bogol- 
ubov quasi-particles with a sound-like spectrum as well 
as gapped spin-one spin- wave excitations. 

To facilitate discussions on the case of Nq = 2k + 1 ± e, 
we introduce, besides spin-one "neutral" excitations, 
"charged" spinless domain- wall excitations. A "charged" 
excitation with a positive or negative charge (±) is de- 
fined by creation operator £>^± and 
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C+ m = ^,- m B+ + = < m B+_ (8) 

where ^v.mW^.m) is a creation (annihilation) operator 
of a spin-one atom with S z = m. Following Eqs.4,8, 
as it is added to one-dimensional DVBCs, a spin-one 
boson fractionalizes into a spin-zero domain wall (with 
an extra atom thus "charged") and a spin- one "neu- 
tral" domain wall because of a two-fold degeneracy at 
No = 2k + 1. This is reminiscent of the physics in con- 
ducting ploymers [20]. At No = 2k + 1 ± e, both types 
of domain walls are free quasi-particles and charged do- 
main walls can condense at the bottom of a band similar 
to that in Eq.7, as a consequence of spin-" charge" sep- 
aration in one-dimensional optical lattices. Spin-charge 
separation was previously pointed out for correlated elec- 
trons in antifcrromagncts [21]. 

Disordered states in square lattices In the limit 
when E CyS are much larger than t, our analysis sug- 
gests that the constrained quantum rotor model (in Eq.l) 
generalized to square lattices should be equivalent to 
an Ising gauge theory in (2+l)D in a quantum spin- 
disordered limit (for detailed discussions on the map- 
ping, see [16]). The effective Hamiltonian is (in an ar- 
bitrary unit) Hp.c. = -r]Tcr x - J2nI\o a z, where a x 
and <j z are Pauli matrices defined at each link and □ 
stands for an elementary plaquette in square lattices in 
this case. The Hilbert space of H p . c . is subject to a lo- 
cal constraint J\ + <j x = (—l) Na at each site, with the 
product carried over all links connected with that sin- 
gle site. In (2+l)D, this Ising gauge field model H p . c . 
is dual to an Ising modcl(frustratcd when No is an odd 
integer) in a transverse field and also represents a parity 
conserved quantum dimer model [22-27]. The frustrated 
Ising gauge field theory has been employed to investigate 
electron bond ordering [24-28] . 

As t becomes much smaller than E s , T approaches in- 
finity. By examining the model straightforwardly, one 
finds that the ground state of H p . c . for A^ = 2k is a non- 
degenerate spin singlet; this is consistent with previous 
discussions on Mott states for an even number of atoms. 
For A*o = 2k + 1 and T ^> 1, the ground state breaks no 
continuous symmetries but does break the lattice trans- 
lational symmetry. It has a four-fold degeneracy and 
supports no gapless modes. This result implies that in 
a square lattice in a quantum spin disordered limit, the 
ground state for an odd number of bosons per site should 
be a valence-bond crystal (columnar-type) state. Two di- 
mensional valence-bond crystals have also been proposed 
in a few recent models on strongly- correlated electrons 
[25-27] and in earlier works [29,30,28,24,31]. 

The ground state \g > on a torus of 2M x 2M sites 
should be topologically identical with one of the following 
configurations 

V G .G.ViM* exp ( 9^e* + (Qo = fa7T.fr))) 10 > • (9) 

»7>£=0,1 



In Eq.9, 5n = 0. .900 = ±5io = 1, S01 = 0, 
h(Q) = exp(iQ • e x ) and g Q0 = ±goi = P 9 10 = 0, 
h(Q) = cxp(«Q • e y ) correspond to four-fold degenerate 
DVBC states in square lattices (^ + (Qo) is defined after 
Eq.4). 

An examination of Eq.9 indicates that both "charged" 
spinless solitons and spin-one "neutral" solitons similar 
to those introduced in one-dimensional lattices interact 
with each other via a linear confining potential in a two- 
dimensional DVBC. The ground state therefore supports 
excitations of S = 2, 1 bound states of soliton anti-soliton 
pairs unlike in one-dimensional DVBC states. 

We should emphasis here the important differences be- 
tween the physics of Mott insulating states of interact- 
ing spin-one bosons and of Hciscnberg antiferromagnets 
(unfrustrated). First, for interacting spin-one bosons, 
the Hilbert space at each site involves different spin 
states^ = 1,3,5,... for an odd A*o much larger than 
unity) unlike in Heisenberg antiferromagnets; this leads 
to a unique Ising gauge symmetry in the problem of inter- 
acting spin-one bosons [9,11], nematic Mott states in high 
dimensional optical lattices [10] and dimerization in low 
dimensional Mott insulating states. Second, phase tran- 
sitions between nematic and spin singlet Mott insulators 
for spin-one bosons can be conveniently investigated by 
changing exchange interactions between two sites or prac- 
tically by tuning laser intensities in optical lattices. This 
opportunity appears to be absent in simple unfrustrated 
Heisenberg anti-ferromagnets. 

Finally, we notice that two-fold degenerate states in 
one dimensional optical lattices could be qubits of a quan- 
tum computer. I would like to thank F. D.M Haldane 
for many stimulating discussions on the BLBQ model, E. 
Demler for an early collaboration on a related subject and 
the year 2001 workshop on "beyond BEC" at Harvard- 
MIT center for ultracold atoms for its hospitality. I am 
also grateful to Bianca Y. L. Luo for encouragement and 
support. This work is supported by the foundation FOM, 
in the Netherlands under contract 00CCSPP10, 02SIC25 
and by NWO Projectruimte 00PR1929. 
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